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Abstract 

The differential cross section of the chargino-neutralino production, qq — > 
X^X } followed by their decays into scalar tau leptons, x^X ~> (T =t i/)(f =F r =t ) 
— > {^^{f^l^vu), is calculated including the effect of spin correlations. In the 
case where f is long-lived, this final state can be fully reconstructed in a hadron- 
collider experiment up to a discrete two-fold ambiguity. Distributions of various 
kinematic variables can thus be observable and tell us about masses and spins of 
superparticles and also parity/CP violation in interactions by comparing with the 
cross-section formula. Observing non-trivial distributions derived in this paper 
will be a good test of supersymmetry. 



1 Introduction 



It is often stated that the LHC is a machine for discovery of new physics and we will need a 
new lepton collider to find out what the actual underlying theory is. This is because most 
of new physics signals at the LHC involve multiple jets in final states which are not simple 
objects to deal with. It is also true that studies of events with missing momentum at hadron 
colliders are more challenging compared to those at lepton colliders because we cannot use 
the momentum conservation in the beam direction. Moreover, unfixed energies of the initial 
partons are another obstacle in studying the exclusive processes. For this reason, most studies 
are limited to forming Lorentz (or boost) invariant quantities out of visible objects to look 
for peaks, endpoints or excesses above expected backgrounds. Such kinds of observables do 
not usually give enough information to determine the Lagrangian parameters. 

Although lepton colliders generally offer a better environment for the studies of exclusive 
processes, at hadron colliders it is not impossible to carry out a detailed study of new-physics 
events if the final states are clean enough. In fact, one of the best-motivated models of new 
physics, supersymmetry (SUSY), may provide such an opportunity. In the case where the 
scalar tau lepton (f) is lighter than the neutralinos and sufficiently long-lived, final states 
of SUSY events have two charged tracks of f rather than a missing momentum associated 
with escaping neutralinos. The presence of such a long-lived charged particle significantly 
improves the capability of the LHC to study SUSY models. 

Although the light f scenario has been treated as an alternative and exotic possibility, it is 
actually neither theoretically exotic nor cosmologically problematic. Since the right-handed 
f carries only the U(l)y quantum number, quantum corrections to its mass through gauge 
interactions are small whereas colored and SU(2) charged sfermions obtain large positive 
contributions. In addition, the Yukawa interaction tends to give a negative contribution 
to the mass. Therefore, it is pretty reasonable to assume that the f is the lightest among 
the superpartners of the Standard Model fields. In such a case, the lifetime of f can be 
very long although the estimate depends on the detail of the model; it can decay into a 
gravitino and a tau lepton through a suppressed interaction if it is kinematically allowed or 
into two Standard Model fermions if i?-parity is violated. There are cosmological constraints 
on such a long-lived charged particle [jjj (H for related works), but those can be evaded as 
long as we do not assume an extremely long lifetime. (See [|| for a recent realistic scenario 
of supersymmetry which predicts a long-lived f and naturally explains dark matter of the 
Universe by gravitinos.) 
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There have been studies of the long-lived f at the LHC, and dramatic differences from 
the stable neutralino scenario have been reported. In Ref. [Q], a technique to reconstruct 
neutralino masses has been proposed by looking for the decay process x° ~^ tt. (See [||, ||] 
for recent studies based on different SUSY models.) A detailed study of measuring the 
mass and the momentum of f in the muon system of the ATLAS detector has been done 
in Ref. || 0, |8[ [|, and it was reported that the mass can be measured with an accuracy of 
0(0.01 — 0.1%) ||. An amusing possibility to collect f's by placing a material outside the 



detectors and measure its lifetime has been proposed in Refs. [10, 11, 12|. Recently, it was 
pointed out that the spin of f can be measured by looking at the angular distribution of the 



pair-production process of f [13|. To discover the long-lived f scenario at hadron colliders, 
various signatures have been considered such as highly ionizing tracks [p~3| , [l~5| , 16], events 
with multiple leptons [|l5|, |16| , and an excess in the dimuon-like events fllq] . (See also [17] for 
a list of various final states.) The usefulness of a px cut (px distribution) in distinguishing a 
f track from a muon has been pointed out in Ref. [18]. 

In this paper, we study the production process of neutralinos and charginos followed by 
their decays into f's. We assume the lifetime of f is sufficiently long (3> ns) so that most 
of the produced f's reach the muon system where their three-momentum can be measured. 
Combined with mass measurements Q, one can reconstruct the four-momentum of the f's. 
We mainly focus on the chargino-neutralino production process since it has the largest cross 
section among the electroweak production processes and the final state is rather simple but 
rich enough to be reconstructed on an event-by-event basis. A particularly interesting process 
is qq — > x^X ~~ > (f ± ^)(f T T ± ) — > (f ± i / )(f ^l^vv), where it is required that the neutralino 
decays into f with the opposite charge to the one from the chargino to avoid a combinatorial 
background. The leptonically decaying r's are selected so that we can easily measure the 
charge of r. The leptonic mode is also cleaner than r-jets with which we need to worry 
about uncertainties such as fake jets and the energy scale. The final state (two opposite-sign 
f's, a lepton and a missing momentum) is clean enough to be compared directly with the 
theoretical calculation. We present a formula of the cross section taking into account the spin 
correlations and demonstrate that various distributions can be seen at the LHC experiments. 
These distributions will be non-trivial tests of SUSY. Methods to measure the neutralino and 
chargino masses by using exclusive processes are also presented. 
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2 Interaction Lagrangian 



There are two types of Feynman diagrams for the x ± X°"P r °duction process. One is through 
an s-channel VF-boson exchange and the others are the t- and it-channel squark-exchange 



diagrams [19]. The interaction Lagrangian for the former diagram is 

C w = xVK-Pl + w R P R ) X -W+ + h.c, (1) 

where wl and wr are coupling constants. We will discuss their relation to the fundamental 
parameters later. For the squark-exchange diagrams, the interaction terms are 

C N = nf WP L u)u{ + nf WP L d)d{ + h.c. , (2) 



Cc = c^{ X +P L u)d\ + cf (x-PLd)u r L + h.c, (3) 

where n^' d ^ and c^*'^ are coupling constants. If we neglect the u- and e?-quark masses, there 
is no chargino coupling to the right-handed quarks. Therefore, only the left-handed (s)quarks 
participate in the diagrams. 

The charginos decay through a term: 

= c^P L v T r^ + h.c. (4) 
There are two terms for the neutralino decay: 

= X°(4 )p * + n ( pP L )rr t + h.c. (5) 

3 The cross-section formula 

We calculate the differential cross section of qq — > x^X° ~^ {f^v)^^^) — > (t v)(f^ I vv) 
in terms of the kinematic variables defined in Fig. [l|. We require the neutralino to decay 
into f with the opposite charge to the chargino, and the r to decay leptonically. From this 
condition and looking at the charges of the lepton and those of two f's in the final state, we 
can tell which f is from the chargino. The angle (0 < 9 < it) is defined as the polar angle 
of the chargino momentum in the center-of-mass (CM) frame where we take the production 
plane to be the x-z plane and the direction of the x-axis is chosen so that the ^-component 
of the chargino momentum is positive. The direction of the z-axis is taken to be that of the 
q momentum. We also introduce angles 6\ and (62 and ^2) which are polar coordinates 



4 




Figure 1: The coordinate systems. 



of the f^ 1 (fT) momentum in the rest frame of (x°) (0 — @i,2 < tt and < (^1,2 < 2-7r). 
Momenta in those frames are related by the following Lorentz transformations: 
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The boost factors are defined by 
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l- x\ + x 2 B I 1 

-7; , PB = \ L 

2x B V Tb 



1b = tt^-Z, (3b = J1~^, (9) 



where xa = my+fy/l and xb = m x o/V§ with s = (P x + + P x ) 2 - F° r later use, we define 

z A = 2xa7a, z B = 2x B jB, (10) 

which are energies of the chargino and the neutralino in the CM frame normalized by yfi/2. 
Finally, we define 

*< = §^, (0<*l< 1), (11) 

where the energies of the lepton (Ei) and r (E T ) can be measured in any frame in the 
approximation m T <C m y o. In this limit, the lepton momentum is pointing in the same 
direction to that of the parent r. 

The cross-section formula can be written in terms of a product of density matrices of the 
production part p ab and the decay parts D A (chargino) and D b B (neutralino) (a, b = 0, 3). 
(See |20| ] for example for methods to calculate the cross section.) By using the narrow width 
approximation, it is given by 

dcos 9 dQi df?2 1 1 9% z aPa ( 1 x ~ 
da = — — dzi ' ~ ' 1 



2 4vr 4vr 1 N c 16vr 2 s \ 1 - 

xB{ X ± -> f ± z^) J B(x° -► f^)^^ -► ^w) 
3 

x £ DXie^p^eMfafazt), (12) 

a,£>=0 

with 

^(0 2 , 02, zi) = ~(1 - «j) [(5 + 5z z - 4z?)D b B (6 2 , fa) - ajv (l + z z - Sz, 2 )^ ] , (13) 

where 52 is the coupling constant of the SU(2)l gauge interaction, and % = mvy/v 7 ! with 
mjy equal to the P^-boson mass. The delta factor is simply 5 b0 = (1,0,0,0). A real-number 
parameter a^v (—1 < cin < 1) represents parity violation in the ^-t-r interaction: 

|„M|2 _ |„( r )|2 

„ _ l"x I K°R I (-i A \ 

N = I (r)i 2 , I (r), 2 - ( M ) 
l n L I + \ n R I 

Once we integrate over the lepton-energy fraction, zi, the Da ■ p ■ Db part reduces to 

Da-P-D b ^Da-P-D b . (15) 
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Note that the term which is proportional to ajv in Eq. fll3| ) vanishes after the integration over 
zi- 

The decay parts D a A and D b B have a simple form: 
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(16) 



where ac is the parity-violation factor in the chargino decay. It always takes the maximum 
value: 



ac = 1, 



(17) 



due to the fact that the neutrinos have only the left-handed chirality. Each component of 
Da and Db corresponds to the expansion coefficient of the Hermitian 2x2 spin-density 
matrices in terms of the unit (a, b = 0) and the Pauli (a, b = 1,...,3) matrices. The non- 
trivial dependencies on angles appear if there is parity violation in the decay vertices. An 
integration over a solid angle df^i (df^) leads to 



D A ■ p ■ D B - p ob ■ D B , (D a - p-D B ^ D a - p a0 D° B ). 
When we perform a further integration of angles, dVtidQ,2, and zi, we obtain 

Da - P ■ D B — >■ p D B — >■ p . 



(18) 



(19) 



The production part p ab is expressed in terms of s, the angle 9 and effective coupling 
factors wl and wr defined by 



w L = w L 



1 1 _ r 2 -(")„(<*)* 
2x1 -V s 92/ V2 
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and 



wr = wr + 



1 1 _ r 2 c W* n W 
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(21) 
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The masses and mg L are those of the left-handed squarks, ul and d B , respectively. The 
components p ab are given by 

P 00 = \(\w L \ 2 + \w R \ 2 )z A z B (l + (3 A f3 B cos 2 9) 
+2Re(w* L w R )x A XB 

T-(\w L \ 2 -\w r \ 2 )zaPacos9, (24) 

P m = 1;(\wl\ 2 + \w R \ 2 )z A x B sin9 
+Re[w* L WR\XAZ B sin 9 

T^(|^l| 2 - \wr\ 2 )z a x b (3 A cos 9 sine, (25) 



p 02 = lm[wlw R }x A z B p B sin9, (26) 



p m = j(\w L \ 2 + \w R \ 2 )z A z B {l + PaPb) cos 9 
+2Re[w* L w r]x ax B cos 9 

t\(\w l \ 2 - \wR\ 2 )z A z B (f3 B + (3 A cos 2 9), (27) 



P 10 = ^(\w L \ 2 + \w R \ 2 )x A z B sin9 
+Re[w* L WR]zAX B sin 9 

T^(|^l| 2 - \wr\ 2 )xazaP a cos 9 sin 9, (28) 



P U = (\wl\ 2 + \wr\ 2 )xax b sin 2 9 

+^Re[w* L w R }z A z B sin 2 9, (29) 



p 12 = --Im.[w* L WR]z A z B p B sin 2 9, (30) 



P 13 = ^(\w L \ 2 + \w R \ 2 )x A z B cos 9 sin9 
+ Re [w* L wr\ za x b cos 9 sin 9 

Tt:(|^l| 2 - \w R \ 2 ) x A z B (3 B sin 9, (31) 



p = lm[w L WR]x B ZAf3 A sin 9, (32) 



p 21 = -lm[w* L w R ]zAZB Pa sin 2 9, (33) 

p 22 = ^Re[w* L WR]zAZB(3Af3B sin 2 9, (34) 

p 23 = — lm[w* L WR]z A XB Pa cos 9 sin9, (35) 

P 30 = ~(\w L \ 2 + \w r \ 2 )z a zb{1 + PaPb) cos 9 
— 2Re[w* L w r]x ax b cos 9 

±- A i\w L \ 2 ~ \wr\ 2 )z a z b (Pa + Pb cos 2 9), (36) 

P 31 = -^(|^l| 2 + \w r \ 2 )z a xb cos 8 sin 9 
— Re[w* L WR]xAZB cos 9 sin 9 

±7}(\wl\ 2 ~ \wr\ 2 )x b z a Pa sin 9, (37) 

p 32 = lnx[w L WR]zBXAfiB cos 9sin9, (38) 

P 33 = ~\{\w L \ 2 + \w r \ 2 )z a z b {PaPb + cos 2 8) 
— 2Ke[w L WR]x A XB cos 2 9 

±\{\w L \ 2 -\w R \ 2 )z A p A cos9. (39) 

With a fixed s, the energies z A (= 2EL±/vI) and zr(= 2E x o/V§) and the velocities Pa and 
Pb are constants as defined in Eqs. (|8|), (||) and (|TQ|) . The spin summed part p 00 has also 



been calculated in Ref. [19] 



By using this cross-section formula we will be able to extract various information such as 
parity and CP violating parameters in the interaction Lagrangian. 



4 Asymmetries vs parameters in the Lagrangian 

The cross-section formula derived in the previous section further simplifies in the case where 
the left-handed squarks are much heavier than the x°/x^ z in the intermediate state, or one 
of the X°/x is Higgsino-like. In such cases, the diagrams with squark exchanges are not 
important and the angular dependencies in wl and wr vanish. In this approximation, 

wl =wl, WR = wr. (40) 
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Figure 2: The p parameter dependence of aw and ^w- The labels A = 1,2 and B = 1 — 4 
represent each mass eigenstate of the charginos and the neutralinos. 



This situation is not unrealistic since quantum corrections tend to make the squarks much 
heavier than other superparticles. In the following discussion we will use this simplification. 
For a more general analysis, one should use the full formula derived in the previous section. 

We have defined two parity-asymmetry parameters ajy and ac(= 1) for the decay pro- 
cesses in the previous section. For the production part p, we define the following three 
quantities: 



aw 



\wl\ 



\wr\ 



I 12 i I 12 ' 
\WL\ + \ W R\ 



2Re[w* L w R } 



Vw 



2lm[w* L w R ] 
\wl\ 2 + \wr\ 2 ' 



(41) 



\wl\ + \ w m 

The matrix p can be expressed in terms of the three quantities and an angle 9. In this section, 
we discuss model parameters and their relations to the observables (ojv, o,Wi £,w-> Vw) defined 
here. 

In the minimal supersymmetric standard model, there are five model parameters which are 
relevant for the process: the Higgsino mass parameter p, the ratio of the vacuum expectation 
values of the Higgs fields tan/3(= (H2) / (Hi)), the gaugino mass parameters Mi and M2, and 
the mixing parameter of the scalar tau leptons 6f (fi = cos Offn + sinflffz,). 

The coupling constants wl,r and n^ R in Eqs. (|l|) and (||) are expressed in terms of mixing 
matrices of the neutralinos Oat, of the charginos Ol and Or, and the mixing angle of the 
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Figure 3: The phase dependence of the CP asymmetry rj\y. The right figure is the same as 
the left figure with a different scale. 



scalar tau leptons Of. The matrices are defined by 



O n M x oO% 



M d j ag 



R M y+ 0\ = M diag 



(42) 
(43) 



where the right-hand-side of the equations are diagonal matrices with real and positive 
eigenvalues. The mass matrices M x o and M x + are 
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The vacuum expectation value u is w = y (.Hi) 2 + (-H2) 2 = 174 GeV. 
In terms of the mixing matrices, the coupling constants are given by 
w L = 92(0* n )b2(01) a1 + ^=(0* n ) B3 (01) A2 , 

WR = 92{0 N )b2{0* R )a1 - ^=(0 N )b4(O r )a2, 



(44) 



(45) 



(46) 
(47) 
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and 



n 



n 



R 



L 




v cos f3 



(O n )b3 cos Of 



(49) 



(48) 



where gy is the coupling constant of the U(l)y gauge interaction. The subscripts for the 
mixing matrices indicate their corresponding components. The indices A(= 1,2) and B{= 
1,...,4) represent mass eigenstates of the charginos and the neutralinos, respectively. The 
second indices of Ol, Or and On are the ones for the interaction eigenbasis; (Wino, Higgsino) 
for charginos and (Bino, Wino, down-type Higgsino, up-type Higgsino) for neutralinos. 

We show in Fig. || the aw and £w factors in a limited case where we fix (tan/3, Mi, 
M 2 ) to be (30, 132 GeV, 250 GeV) and vary the fi parameter from 150 GeV to 350 GeV. 
The labels A = 1, 2 and B = 1, ...,4 represent each mass eigenstate of the charginos and the 
neutralinos, respectively. For a small value of /i, the lighter chargino (A = 1) is Higgsino-like. 
As \jl increases, the lighter chargino goes through the mixed region (/i ~ M 2 ) to the Wino-like 
region (/i 3> M 2 ). The largest cross section is for A = 1 and B = 2 which gives aw — 
and £w — 1- in the two extreme limits where the produced chargino and neutralino are both 
purely Higgsinos or both purely Winos, there is no parity violation in the interaction vertex 
(i.e., aw = 0) since they are vector- like particles. A large deviation from aw ~ is possible 
when there is a significant mixing among the Higgsinos, Wino and Bino. 

The CP asymmetry r]w is calculated with varying the phase of the \i parameter in Fig. ^. 
The parameters are fixed as (tan/3, Ml, M 2 , \fi\) = (30, 132 GeV, 250 GeV, 200 GeV). A 
large CP asymmetry is obtained only for A = 2 and B = 1. Other asymmetries are at most 
of order a few percent. (Note that the relative phase between M\ and M 2 , arg(MiM|), is 
also an independent physical parameter.) 

Finally, the parity asymmetry oat in the neutralino decay is calculated with (tan/3, Mi, 
Ms) = (30, 132 GeV, 250 GeV) in Fig. [|. We take three values of the f-mixing parameter 9? = 
0,7r/4 3 vr/2. Fig. | shows that oat is highly dependent on the model parameters; specifically 
the properties of the neutralino and f. This parameter does not necessarily vanish in the 
pure Higgsino or Wino limits due to the chiral nature of the tau lepton. For example, if the 
neutralino is purely Higgsino and the stau is left-handed (right-handed), the tau lepton must 
be right-handed (left-handed), i.e., a^r = — 1 (+1) ; although we need to take into account 
mixings in more realistic cases. 

We can see that the asymmetry parameters, especially a^, are sensitive to the model 
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Figure 4: The \x parameter dependence of the parity asymmetry am. 

parameters. As we will see in the next section, the parameters, a at, aw, VW can be measured 
by looking at asymmetries in various distributions. Measurements of those asymmetries 
together with the mass measurements provide us with information on combinations of the 
neutralino/chargino mixings and the f mixing. 

5 Angular and energy distributions 

We can obtain various one-dimensional distributions by integrating over the remaining vari- 



ables in the differential cross section (12). Here we adopt the simplification that the squark 
diagrams are not important. 

We list in Table [l] transformation properties of the angles under the charge conjugation 
(C), the parity transformation (P), and the time reversal (T). The transformation properties 
of asymmetry parameters are assigned in the right table. With these assignments, the 
interaction Lagrangian in Eqs. (Jll — 1ST) are "formally" C, P, and T invariant. These are helpful 



in understanding the resulting distributions. 
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Table 1: C, P, and T transformation properties of angles denned in Fig. || and of the 
asymmetry parameters. With these assignments in the right table the interaction Lagrangian 
in Eqs. (|l]-[|) are "formally" C, P, and T invariant. The symbol ± represents the charge of 
the chargino. 



Zi distribution 

By integrating over 9 
da 



n, 91 



), and 



), we obtain 



a(qq X ± X°)B(x ± -> r ± u)B( X ° -> f^r^B^ -» l ± vv)dz l 
X-(l - Zl ) [(5 + 5zi - 4zf) - a N (l + z x - 8zf)] . 



(50) 



This is the well-known polarization dependence of the lepton-energy distribution occurring 
in leptonic r decays pj| . Since z\ is a rotation and boost invariant quantity (in the limit of 
m T /m x o <C 1), we can measure this distribution in the laboratory frame. This distribution 
will tell us about the parity asymmetry a^ in the neutralino decay through the r polarization. 
This distribution will remain unchanged when we include the squark diagrams. 



cos Ox distribution 

By integrating over 9, 4>i, (92, 92), and zi, we obtain 

da = a ( qq - ^ x ± x °)B(x ± - f ± v)B( X Q - f T r ± ) J B(r ± -> \±vv) 
x [1 + awacfi(/3A, /3b) cos 9i] . 

Recall that ac = 1. The function /i is given by 

3/3a + /9s 



d cos 9\ 



3 + fofo + 3^ x /(l - )(i - 



(51) 



(52) 



Because this cos 9\ distribution is P-even (see Table |IJ) , a non-trivial distribution requires 
parity violation in both the production process (aw) an d in the decay of the chargino (ac)- 
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In hadron collisions, event rates are obtained after a convolution with the parton distri- 
bution functions. The actual distribution is oc 1 + aw (fi) cos 0i where (f\) is an averaged 
value of fi. The value of f\ vanishes in the threshold production limit ((3a —> 0, (3b — ► 0), 
and it approaches unity for a boosted event ((3a — ► 1, (3b !)■ (A larger asymmetry 
can be observed if we select events with large s, although the number of events decreases 
exponentially if the lower cut on s is increased.) Observing this asymmetry will provide 
evidence of both the chargino spin and of parity violation in the weak interaction of the 
charginos and neutralinos. 

cos 02 distribution 

A similar distribution is obtained when we integrate 9, (9±, 0i), 02, and zf. 

da = a( q q^x ± X°)B(x ± ^^u)B( X ^^^)B(r ± ^l^) d -^ 

x [1 + a w a N f 1 ((3 B , (3a) cos 6 2 ]. (53) 

The angular dependence is due to both the neutralino spin and the parity violation occurring 
both in the production and the decay of the neutralino. 

cos 0! cos 6 2 distribution 

A non-trivial correlation is present between the angles on both sides of decays. The 6\ and 
O2 dependence of the cross section is 

da = a(qq^x X )B(x r v)B( X -► r+r )B(t -► I vv) 

x [1 + a w acfi((3A,(3B) cos 0i + a w a N fi((3 B , (3a) cos 9 2 

+aca,N f2((3A, (3b) cos6»i cos0 2 ], (54) 



where 



3(3aPb + 1 + 

I2((3a,Pb) = V , (55) 

3 + (3 A (3 B + ^w^(l-l3 A )(l-(3 2 B ) 



Integrating over 9\ and 02 keeping the product cos0i cos02(= y) fixed, we obtain 
da = a(qq - xVMx* - r ± ^)5( X ° - f ^r ± ) J B(r ± - Z±^) 



2 

x [1 + a c a N f 2 (PA,PB)y]log\y\. (56) 
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The non-trivial part (the second term) is due to the spin correlations between the chargino 
and the neutralino. Parity violation (a^v 7^ 0) biases the distribution towards a positive or 
negative value of y. This distribution is independent of the asymmetry parameter aw m 
the production process. Note also that the function fa does not vanish in the limit of the 
threshold production (/2 — ► 1/3), although it is maximized in the boost limit (/2 — > 1). 
Confirming this correlation will be an interesting test of the model. 

</>! distribution 

A non-trivial distribution of the azimuthal angle <f>x takes place due to parity violation in the 
chargino decay: 

da = a{qq^x ± X Q )B{x ± ^T ± v)B( X ^fT T ± )B{T ±^ l ± v „ ) ^± 



x 



7T 2 ,„ _ vr 2 



1 ± T^ a c9i(PA,PB) cos^i ± — ocfJ7w52(/3 J 4)^B)sin<^i 
lb lb 



(57) 



where 



9i(Pa,Pb) = * ; V (58) 



1 + p A p B /3 + i w J(l - (3 2 A )(1 - (5 2 B 



PaJi-P 2 b 

92Wa,Pb) = V , =■ (59) 

1 + p A p B /3 + t w ^J(l - (3\){l - f3%) 

The (pi dependence appears even if aw = 0. This is somewhat surprising once we realize 
the fact that cos^i is P-even and the chargino decay violates parity (ac 7^ 0). In order for 
the distribution to be formally P-invariant there should be another interaction that violates 
parity. This is in fact supplied by maximal parity violation in the weak interaction of the 
quarks in the production process. This fact means that to observe the distribution one needs 
to measure the direction of the quark (or the anti-quark) . This conclusion can be also seen in 
Fig. ffl, because knowledge of the quark direction is necessary to define the angles <p\ and <fo- 
The different signs for the x + X° an d X~X° productions can be understood by the fact that 
cos 0i and sin^i are CPT-odd. The sin^i dependence (phase of the (pi oscillation) measures 
CP (or T) violation in the production process, r\w ■ 

The function g\ is maximized at the threshold limit and vanishes in the boost limit, in 
contrast to the case of the polar-angle dependencies. In the threshold limit, the coefficient of 
cos (f)\ is 7r 2 /lb. The CP asymmetry vanishes in both the threshold and the boost limits. 
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The determination of the spin of the intermediate particles by looking at the azimuthal- 
angle distributions (frequencies of the cj> oscillations) has been discussed recently in Ref. [f2^]f|. 



4>2 distribution 

A similar distribution is obtained when a at is non- vanishing: 



da = a(qq - X ± X°)B( X ± - f^Btf - f^r^B^ - /±z^, 

1 T T7:awgi(/3B, A-l) cos 02 =F TF a NVwg2(PB,PA)sm^2 
lb lb 



2 



(60) 



Other distributions 



Although we will not study them in this paper, there are various kinds of other non-trivial 
distributions. For example, the distribution of the difference of the angles (j>\— 4>2 also depends 
on the CP-violation parameter rjw- In the reconstruction of this angle at hadron colliders we 
do not need to know the direction of the q or q in the initial state in contrast to the case of 
the angles 0i and fo- This is an advantage especially at pp colliders. Analytic formulae of 
such distributions can easily be obtained from the full cross-section formula. 



6 LHC studies of x + X an d X^X productions 

In this section, we demonstrate a possible strategy for the study of the production processes 
of charginos and neutralinos by performing a Monte Carlo simulation. We use the following 
simplified model for generating events: 

H = 300 GeV, Mi = M 2 = rriu L = m~ d = 5000 GeV, (bl) 



tan p = 10, m fL = 5000 GeV, m fR = 100 GeV. (62) 

With this choice of parameters, all the SUSY particles decouple from low energy except for 

the Higgsinos and the right-handed f. The chargino and two light neutralinos are purely 

*As we have seen above in the case of the X ± X° production, parity violation is needed at both the production 
and the decay vertices in order to develop a (j>\ or 4>2 azimuthal-angle dependence. These conditions are in fact 
general requirements for 2-^2 fermion pair production with subsequent two-body decays of each fermion. 
Therefore, the method of Ref. [^] should work only in a limited case. For example, there is no azimuthal-angle 
dependence of the differential cross section in processes where fermions are pair produced through QED or 
QCD interactions such as tt pair production or the production of a gluino pair at hadron colliders (unless the 
beam is polarized). In such processes, the angular correlations between two decays, such as the distribution 
of cos 9i cos 62 or cj>i ±02, will instead be useful if there is parity violation at the decay vertices. 
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Higgsino-like and the masses are calculated to be: 



m fl = 109 GeV, m x + = 300 GeV, m x o = 299 GeV, m x o = 301 GeV. (63) 

Although there are two mass eigenstates for the neutralinos due to a small mixing with 
gauginos, they almost behave like a single Dirac fermion. We do not distinguish x± an d X2 
in the following analysis. The lifetime of f\ is assumed to be much longer than the typical 
collider time scale (l/IV 3> ns). The branching fractions of the chargino and the neutralino 
decays are of course, 

B{x ± -> ^v) ~ 1.0, B(x° -> fV 1 ) ~ 0.5. (64) 

The asymmetry parameters in this model are calculated to be 

a N = 1.00, a w = 0.00, £ w = 1.00, r) W = 0.00. (65) 

Note that this parameter choice is simply for a demonstration and not particularly motivated 
by any fundamental model which realizes a light f. We use this model as the first trial of 
the study of production events of charginos and neutralinos in the long-lived f scenario. In 
more realistic situations, other mass eigenstates will be produced which contaminates the 
analysis of the leading production process. Since the importance of such effects depends on 
the detailed structure of the models, we use the above clean model as a toy example. The 
values of asymmetry parameters, aw = 0.00 and rjw = 0.00, are not very interesting ones, but 
in fact, as we see later we need to first study these trivial cases in order to confirm whether 
there is a fake distribution caused by false solutions, which appear in the reconstruction of 
kinematic variables. In order to measure the asymmetry, we need to understand whether a 
non-trivial distribution is fake or physical. 

We have generated 26,000 events of the electroweak production processes of SUSY parti- 
cles (including the f pair-production process) in the pp collision at = 14 TeV by using the 
Herwig 6.5 event generator |f23[|. The spin correlations have been implemented for the x + X° 



production and their decays [24, [25|]. This number of events corresponds to an integrated 
luminosity of 100 fbr 1 at the LHC. (We will use 300 fb _1 of data for some of the analysis 
of angular distributions.) We have used the CTEQ5L library |26j for the parton distribution 
function. For the r decay, we have used TAUOLA 2.7 package ]27j so that the spin information 
is maintained. A detector simulator AcerDET 1.0 [^q] has been used for the event analysis. 

In the following analysis, we assume that the mass of f is known by the method of Ref. ||, 
and we ignore the resolution of the f-momentum measurements which is of order a few percent 
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Figure 5: The Mti distribution of the chargino-pair production. The chargino mass can be 
extracted by looking at the endpoint. 

in the ATLAS experiment j?]]. One should note that the accuracies of the measurement quoted 
below are somewhat optimistic for this reason. We also assume perfect efficiencies of the f 
identification and of the f-charge measurement for f tracks with px > 10 GeV and \rj\ < 2.5. 
By requiring two f's, there is no Standard Model background with this assumption although 
in actual experiments one needs to take into account mis-identifications of muons as f. 

We first discuss possible methods to measure the masses of the chargino and the neutralino 
through the exclusive production processes. Measurements of the asymmetries from looking 
at the angular and energy distributions studied in Section [| will then be demonstrated. 

6.1 Chargino mass determination by chargino-pair production 

We present a method to measure the chargino mass exclusively from the chargino pair- 
production process. The final state of the process is two opposite-sign f's and missing 
momentum from the two neutrinos. 

Although we cannot reconstruct the chargino momentum on an event-by-event basis, the 
endpoint analysis developed in Ref. [ |29| | can be used to extract the chargino mass. The 
method is to form a quantity Mt2 defined by 

M? 2 = min [max{m^(p Tf - , PivJ, r?4(p rf +, ptv 2 )}] , (66) 

Pt i , 1 +PT,^ 2 =Pt 1SS 
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where Pxf- and PTf+ are the transverse momentum of the two r's and p™ lss is the missing 
transverse momentum. The transverse mass txit is defined by 

m\ = {Et? + E T uf ~ \PTf + PTu\ 2 

= m 2 . + 2(E Tf E Tl/ -p T fp Tu ), (67) 

where 

E\ = m 2 + |p T | 2 . (68) 

The quantity Mt2 is designed to have the endpoint at the mass of the intermediate particle. 

In order to select the x + X~ events, we have imposed the following jet and lepton vetoes: 

Nf+=N f - = l, Njipr > 30 GeV) = 0, Ni(p T > 6 GeV) = 0. (69) 

We do not need to impose a tight cut on the missing momentum since f + f~ events do not 
contribute near the endpoint of the distribution. The Mt2 distribution is shown in Fig. ||. 
There is a clear endpoint around the input chargino mass, 300 GeV. By fitting with a linear 
function, we obtain the endpoint: 303.2 ± 0.7 GeV, which is slightly larger than the input 
value due to the resolution of the missing transverse momentum]]. 

6.2 Neutralino and chargino mass determination by chargino-neutralino 
production 

We show in this subsection that quite accurate measurements of the neutralino and chargino 
masses are possible by analyzing exclusive processes. Combining various methods described 
below, we will be able to measure the masses at the level of a few GeV. 

6.2.1 Endpoint analysis for the neutralino mass 

The neutralino mass can be measured by looking for an endpoint of the invariant mass 
distribution of the f^l^ pair from the neutralino decay followed by the leptonic tau decay. 
The x x°-production process can be selected by requiring two f's and an isolated lepton: 

N f +=N f -=l, Njipr > 30 GeV) = 0, N t (p T > 10 GeV) = 1. (70) 

As we discussed before, we require that two f's have opposite signs so that there is no 
ambiguity in selecting f from the neutralino decay. 

The invariant mass distribution of the f^l^ pair is shown in Fig. ^. An accurate 
measurement of the neutralino mass is possible by this method (300 ± 3 GeV). 
^One should use a Gaussian-smeared line to take into account the effect of finite resolutions. 
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Figure 6: Invariant mass distribution of f =F Z =t: pairs. The endpoint shows the neutralino mass. 
6.2.2 Solvability analysis for the neutralino mass 



By using the information of the chargino mass measured by the X + X P a h production process, 



we can obtain the neutralino mass by a similar method proposed in Refs. [30, |3l]]. (See also 
[32] for a similar analysis for the measurement of the top-quark mass in the di-lepton events 
from the tt productions at the LHC.) Since the final state is relatively simple, we can solve 
the kinematics on an event-by-event basis by postulating a neutralino mass. By maximizing 
the solvability (number of events which can give physical solutions normalized by the total 
number of events analyzed), we can obtain the correct neutralino mass. 



The equations to be satisfied are 

(P f ± + P v f 



(71) 



zi x 



px _|_ \ px px 

r v ' r l± -"miss' 



py _i_ _ z ± pv _ pv 

1 v ~ 1 l± 1 miss' 



whereas there are four unknowns in these equations: 

px py pz 



(72) 
(73) 
(74) 

(75) 
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Figure 7: The solvability analysis for the neutralino mass. The right figure shows the 
solvability near the peak. 



Equation ([72]) is a linear equation for z\ in the approximation of m\ = 0, and by using 
the solution, Eqs. ( |73[ ) and ( f74] ) become also linear equations for P^ and P;, respectively. 
Eq. ([Tip, on the other hand, is a quadratic equation for P* , and therefore there can be either 
zero or two real- number solutions. (In general, an equation of this type may have a unique 
physical solution by the constraint E v > 0. However, one can show that Eq. (|7l| ) always 
have zero or two solutions by using the fact that E f ± > |Pf ± |.) The number of events should 
be maximized at the correct neutralino mass when we impose conditions: < z\ < 1 and 
existence of real-number solutions of Eq. (71). 

The number of events with a physical solution is shown in Fig. ^ for various input 
neutralino masses. It indeed shows a sharp peak at the correct neutralino mass, 300 GeV. We 
have used the correct value of the chargino mass in the analysis. In the actual situation, the 
experimental error in the chargino mass will propagate into the error in the peak location. 
Fitting with two linear functions near the peak, we find that the neutralino mass can be 
measured quite accurately (301 ± 2 GeV) if the chargino mass is known. 

Note again that this analysis is not completely realistic. We have ignored the momentum 
resolutions of f tracks and assumed the perfect identification efficiency. We leave more 
realistic studies to future work. 
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Figure 8: The transverse mass distribution for the chargino mass measurement. 



6.2.3 Transverse mass analysis for the chargino mass 

Once we know the neutralino mass by, for example, the method of the endpoint of the 
Mj^i± distribution, the transverse momentum of the neutrinos from the chargino decay can 



be reconstructed without the two- fold ambiguity from Eqs. (72-74). We can then form a 
transverse mass in Eq. ( |67| ) and the chargino mass can be obtained by looking for an endpoint 
of the distribution. 

We show in Fig. || the distribution of the transverse mass, My. We can see a sharp 
peak near the correct chargino mass. The endpoint is again smeared by the resolution of the 
missing transverse momentum. An appropriate fitting is necessary for the extraction of the 
chargino mass. For a simple fitting by a linear function, we obtain a significantly larger value 
(309.2 ± 0.8 GeV) due to the finite resolution. 



6.2.4 Solvability analysis for the chargino mass 

The solvability analysis can also be done for the chargino mass once we know the neutralino 
mass. The solvability is plotted in Fig. |9] where we see that the solvability saturates near the 
chargino mass. 

By looking for a point where the solvability saturates, we can obtain the chargino mass 
(303 ± 1 GeV). 
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Figure 9: The solvability analysis for the chargino mass measurement. The right figure is the 
same analysis with a better resolution near the threshold. 

6.3 Energy and angular distributions 

Now we examine whether the energy and angular distributions obtained in Section [| are 
visible in actual experiments. An especial concern is that there is always a false solution 
in the Eqs. (ffl|-74), which may destroy the theoretical distributions. One purpose of this 



analysis is to understand the effect of the false solution. 

In the analysis, we have used the events passed through the selection cut in Eq. (|70|). We 



assume in the following that the chargino and neutralino masses are known and ignore errors 
in the mass measurements. 



6.3.1 Zi distribution 

This distribution measures the polarization of the r lepton from the neutralino decay. We do 
not need to distinguish events with different lepton charges since the theoretical distributions 
are the same (Eq. (f>0"l)). 

The measurement of the energy fraction of the lepton, zi, does not suffer from the two-fold 
ambiguity since Eq. ( |72[ ) is a linear equation in z\ in the approximation of mi = 0, and we do 
not need to know P*. 

The distribution is shown in Fig. pi]. We fit the distributions with the theoretical curve 
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Figure 10: The reconstructed distribution of the lepton-energy fraction (z{) . The solid curve 
is the best fit with a theoretical function in Eq. fl50|). The curves with a at = 0, —1 are also 
shown (dashed curves). 

in Eq. (^) by making a^v a parameter. We obtain a^v = 1-1 ± 0.2 (solid curve) for = 1.0 
in this model. Curves with oat = and on = — 1 are shown in the figure (dashed lines, we 
used the same normalization with the solid curve). We can see that the best-fit curve can be 
discriminated from those models. The region with small z\ is affected by the px cut on the 
leptons. This region is omitted from the fitting. 



6.3.2 cos #i distribution 



This distribution measures the parity asymmetry in the x^X° production process, aw, 



through Eq. (51). The averaged value of the function fi in Eq. (52) weighted by the cross 
section depends on a selection cut on s (s = {P x ± + P x o) 2 ). It is an increasing function of 
s m ; n , but the number of events rapidly decreases with s m i n . In the model we simulated, the 
cross section falls off as 

da(s) 



oc exp 



-4.3 



TeV 



(76) 



The averaged value defined by 



(/l(Viniin)) 



da(s) 



fi(\fl)d\/J 



da(s) 



dVI 



(77) 
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Figure 11: The reconstructed cos#i (left) and cos 62 (right) distributions. The false solutions 
are included. A selection cut \f§ > 900 GeV is imposed on both of the solutions in each 
event. 

is then estimated to be, for example, 

(/1 (600 GeV)) = 0.52, (/i(900 GeV)) = 0.74, (78) 

for m x + ~ m x o ~ 300 GeV and £w = 1- Since the asymmetry will be diluted by false 
solutions as we see below, it is necessary to impose a cut s m i n in order to expect a large 
asymmetry. 

The angle cos 9\ which is defined in the rest frame of the chargino is expressed in terms 
of the f energy in the CM frame: 

E CM _ E (}) 

cos e 1 = -^ — (79) 

where ja and (3a are defined in Eq. (@), and E { }} and pj£ are the energy and momentum in 
the rest frame of the chargino, respectively. They are given by 



eV = ^^, m = \l(m 2 -ml (80) 



2m x+ ' f± VV 

When we calculate the f energy in the CM frame from the quantities measured in the 
laboratory frame by boosting to the z-direction, one encounters the two-fold ambiguity for 
Py m Eq. ©. 
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In order not to develop a fake distribution caused by the false solution, we need to be 
careful when we impose a selection cut on s. We have examined the following three methods 
of imposing a <§min cut. One is to choose a solution which gives a smaller value of s, and 
impose a selection cut vi > 900 GeV on the chosen event. This strategy effectively picks up 
the true solution (with the probability of about 63%) because of the distribution in Eq. (|76l). 
However, this method causes a bias in the cos#i distribution towards larger cos#i. For each 
event, it is likely that the solution with larger cos#i (which would mean that the neutrino is 
emitted to the opposite direction to the chargino in the CM frame) gives a smaller value of 
s, and thus such a solution is more probable to be chosen. This correlation causes bias. 

The next strategy is to use all the solutions with \/s > 900 GeV. That is, if we have two 
solutions which satisfy the cut in an event, we use both solutions. If there is only one solution 
with V§ > 900 GeV, we use that one. This strategy causes a fake distribution towards smaller 
cos#i this time. Since we impose a lower cut on s, this strategy tends to select a solution 
with larger s. By the same reason as above, this tends to pick up a solution with smaller 
cos 6>i . 

The above two lessons lead us to a good strategy to avoid the bias. It is to use both 
solutions in each event and impose an s cut on both of the solutions, i.e., we throw away an 
event if there is a solution with \fl < 900 GeV even though it may be a false solution. By 
doing that, the probability of selecting the true solution is exactly 50%, and there is no obvious 



reason to expect a fake distribution. We show in Fig. 11 the reconstructed cos 6\ distribution 
by using the strategy (left figure). A flat distribution is obtained which is expected in this 
model because aw = 0.0. For a more general case, the slope of this distribution should give 
approximately a^(/i)/2 where the factor of two comes from the effect of false solutions. As 
the statistical error of p2 in Fig. |ll] is ±0.07, the establishment of aw 7^ at the 3a level 
would require |aiy| > 0.6. 

6.3.3 cos #2 distribution 

This measures the product of the parity asymmetries a^ and aw in Eq. (^) through the 
spin correlation of the neutralino. By the same strategy as in the cos 9\ case, we obtain a flat 



distribution for cos #2 as expected. It is shown in the right panel of Fig. 11 
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Figure 12: The reconstructed w(= cos 6\ cos 62 (log | cos #1 cos 62 \ — 1)) distribution. A 
deviation from the flat distribution indicates the presence of the spin correlations and parity 
violation in the chargino and the neutralino decays. The false solutions are included. A 
selection cut \/s > 900 GeV is imposed on both of the solutions in each event. The right 
figure is the same but with 300 fb _1 of data. 

6.3.4 cos #1 cos 62 distribution 

Although the cos 0\ and cos 62 distributions are trivial in this particular model due to a\y = 
0.0, there can be a non-trivial correlation between cos#i and cos 62- An example is the 
distribution of the product cos 0\ cos 62 which measures the product of the parity asymmetries 
in the neutralino and chargino decays independent of aw (see Eq. ([56])). 

We define a variable, 

w = h(y) = y(log \y\ — 1), y = cos 9\ cos 62- (81) 
The theoretical distribution in Eq. ( |56|) in terms of the variable w is 

da oc (1 + a N (f 2 )h- 1 (w))dw, (82) 

where h^ 1 is the inverse function of h(y), i.e., y = h" 1 {w), and —1 < w < 1. The w 
distribution is flat in the parity conserving case (a^ = 0). The deviation from the flat 
distribution is a signature of parity violation. The averaged value of (/2) depends on V§ m i n : 

(/ 2 (600 GeV)) = 0.58, (/ 2 (900 GeV)) = 0.75. (83) 
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Figure 13: The reconstructed <p\ distributions for the X + X° (left) and the \ X° (right) 
production events. 

Therefore, with the strategy for the selection cut discussed before, we expect an asymmetry 
a Af(/2)/2 — 0.38 in this model, where the factor of two is the effect of fake solutions. 

The reconstructed w distribution is shown in Fig. ^ where we see deviation from the flat 
distribution. The right figure is the same analysis with 300 fb _1 of data. We fit the histogram 
with the function in Eq. ( |32"| ) and obtained a significant asymmetry, a_/v(./2)/2 — 0.47 ± 0.12 
(0.49 ± 0.07) which deviates from zero by 4a (7a) with 100 ftr 1 (300 ftr 1 ) of data. A 
somewhat larger value compared to the expectation (0.38) can be understood by the fact 
that the effective vlmiii is larger than 900 GeV because we have imposed a cut on both of 
the solutions. 

Observation of this distribution together with the measurement of a at by the Z\ distribu- 
tion will be a quite interesting confirmation of the spin-spin correlations. 

6.3.5 0i distribution 

Non-trivial azimuthal- angle distributions show up when there is parity and/or CP violation in 
the decay vertex. In order to measure this, we need to completely reconstruct the kinematics 
such as the angle 9. The angle <j)\ is expressed in terms of the angle 9 and the three- momentum 
of in the CM frame: 

Wi = 77^ (P P C m^ (O^'i^)' (8 4 ) 

(-Pf)cM cos 9 - {P~)cm sin 6 
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with 

W = 0i (if (-Ff)cM > 0) 



1 = #+7T (if(P|) C M<0) • (85) 

In order to define the CM frame in Fig. [l], we need to know the direction of the anti-quark 
which can be determined only statistically in pp-collision experiments. We take the z-direction 
to be the same direction as that of the total momentum, P z = P* ± + P* , in the laboratory 
frame since the q parton tends to carry a smaller momentum. In order to reduce the number 
of mis-choices, we impose a cut: P z > 1200 GeV. 



The (pi dependent part of the distribution in Eq. (57) has opposite signs for x + an d x~ 
productions. We do not impose a cut on s because the g± function in Eq. (|5^) takes its 
maximum value at the threshold production. (In order to look for a CP asymmetry, it may 
be better to impose a cut. See Eq. (|59|).) We also use both solutions for P„. The averaged 
value of the functions g\ and g2 are: 

2 2 

^<0i)=O.51, ^( 52 )=0.16. (86) 

We expect that these values will be affected due to the existence of the false solution. 

The distributions are shown in Fig. 13 with 300 fb _1 of data. The left figure is the 



distribution of the x + X° events (i.e., the events with a positive-charge lepton) and the right 
figure is from the x~X° events. We fit the histogram by a function: 

p{<t>i) =Pl(l +P2cos0i +P3sin0i). (87) 

A qualitatively correct behavior is obtained in the x + X° events, i.e., p2 > and p% = 0, but 
X X° events do not show the expected behavior of p2 < and P3 = due to poor statistics 
and the selection cut on P z . We can see from the figures that the selection cut on P z tends 
to give a fake distribution peaked near (pi ~ and 2ir for both x + X° an d X~X° events. One 
may be able to avoid this by imposing the P z cut on both solutions as we have done in the 
study of the cos#i distribution, but it significantly reduces the statistics. A looser cut on 
P z results in a fake distribution by the mis-choice of the z-direction. Nevertheless, it is not 
a problem for observing a non-trivial distribution since the theoretic distribution is different 
for x + X° an d X~X° productions. For example, one can try to rescale the histogram of the 
X"X° events and subtract from (or add to) that of the x + X° events in order to eliminate (or 
understand) the fake distribution. 
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Figure 14: The reconstructed 02 distributions for the X + X° (left) and the \ X° (right) 
production events. 

6.3.6 02 distribution 

The (p2 distribution can be obtained by the same method. The distribution clearly shows an 
expected behavior in Eq. (|60|) with cin = 1.0 and rjw = 0.0. 

7 Summary 

If f is the lightest among the superpartners of the Standard Model particles, the SUSY 
signatures at the LHC experiments will be very different from the stable neutralino scenario. 
We have demonstrated that the production processes of the neutralinos and charginos have 
rich information on model parameters. The spin correlations of intermediate particles give 
rise to interesting non-trivial distributions in various kinematic variables. In previous studies 
of SUSY models at hadron colliders, the production of neutralinos and charginos have been 
usually thought of as good processes to discover SUSY through multi-lepton final states. In 
the stable neutralino scenario it is nevertheless challenging to extract information on models 
out of those processes because of the small cross sections and difficult kinematics due to 
escaping neutralinos. Much attention has been paid to production processes of colored 
superparticles and their cascade decays for the measurements of the model parameters. 
However, as we have shown, the chargino-neutralino production process may provide us with 
the best opportunity for understanding SUSY models in the long-lived f scenario. 
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The study presented in this paper is not fully realistic in several senses. We have not 
included the momentum resolution of the f tracks or efficiency of the identification. Also, in 
the study of various distributions, we have ignored errors in the measurements of the chargino 
and the neutralino masses. We have used the transverse missing momentum evaluated by 
the fast detector simulator, but the resolution may be very different in real experiments. The 
trigger efficiency of the process has also been ignored. A more detailed analysis is necessary 
when we discover the long-lived f . The analytical formulae presented in this paper will be 
useful in such future studies. 



This work is inspired by studies of the electroweak theory in Refs. [33] where the dif- 
ferential cross section of the process e + e~ — * W + W~ is calculated including the effect of 
spin correlations. These various distributions are studied for the purpose of confirming the 
SU(2)l x U(l)y gauge interactions. The density matrix calculated there has been used to put 



constraints on anomalous interactions among gauge bosons at the LEP-II experiments 34]. 
If f is long-lived, the cross-section formula calculated in this paper can be used as a good 
test of SUSY at the LHC experiments just like we have confirmed the Standard Model at the 
LEP experiments. 

We here comment on the x°X° production processes which we did not study in this paper. 
In many cases, these processes have smaller cross sections than the x^X process. Since the 
Z-boson vertex involving the same mass eigenstates, Z — Xi — Xii vanishes identically, the 
main production process is x^Xj with i ^ j. If we require the opposite charges for two f's and 
the leptonic decays for both of the r leptons, the number of events will get smaller. However, 
for the reconstruction of the kinematics, it is much simpler than the x^X production events. 
We can reconstruct the final state without the knowledge of the neutralino masses. Moreover, 
there is no discrete ambiguity for the reconstruction. The study of these processes will also 
be important if f is long-lived, although it may be challenging due to the limited statistics. 

Acknowledgments 

I thank Michael Graesser for reading the manuscript, stimulating discussions and useful 
comments. I also thank Alex Friedland for discussions and useful comments. 



32 



References 

[1] T. Gherghetta, G. F. Giudice and A. Riotto, Phys. Lett. B 446, 28 (1999) [arXiv:hep- 
ph/9808401]; T. Asaka, K. Hamaguchi and K. Suzuki, Phys. Lett. B 490, 136 (2000) 
[arXiv:hep-ph/0005136]; J. L. Feng, S. Su and F. Takayama, Phys. Rev. D 70, 075019 
(2004) [arXiv:hep-ph/0404231]; M. Pospelov, Phys. Rev. Lett. 98, 231301 (2007) 
[arXiv:hep-ph/0605215]; K. Kohri and F. Takayama, Phys. Rev. D 76, 063507 (2007) 
[arXiv:hep-ph/0605243]; M. Kawasaki, K. Kohri and T. Moroi, Phys. Lett. B 649, 436 
(2007) [arXiv:hep-ph/0703122]; J. Pradler and F. D. Steffen, arXiv:0710.2213 [hep-ph]; 
K. Jedamzik, JCAP 0803, 008 (2008) [arXiv:0710.5153 [hep-ph]]; M. Kawasaki, K. Kohri, 
T. Moroi and A. Yotsuyanagi, arXiv:0804.3745 [hep-ph]. 

[2] T. Moroi, H. Murayama and M. Yamaguchi, Phys. Lett. B 303, 289 (1993); E. Holtmann, 
M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 60, 023506 (1999) [arXiv:hep- 
ph/9805405]; J. L. Feng, A. Rajaraman and F. Takayama, Phys. Rev. Lett. 91, 011302 
(2003) [arXiv:hep-ph/0302215]; J. L. Feng, A. Rajaraman and F. Takayama, Phys. Rev. 
D 68, 063504 (2003) [arXiv:hep-ph/0306024]; J. L. Feng, S. f. Su and F. Takayama, 
Phys. Rev. D 70, 063514 (2004) [arXiv:hep-ph/0404198]; M. Kawasaki, K. Kohri and 
T. Moroi, Phys. Lett. B 625, 7 (2005) [arXiv:astro-ph/0402490]; M. Kawasaki, K. Kohri 
and T. Moroi, Phys. Rev. D 71, 083502 (2005) [arXiv:astro-ph/0408426]. 

[3] M. Ibe and R. Kitano, JHEP 0708, 016 (2007) [arXiv:0705.3686 [hep-ph]]. 

[4] I. Hinchliffe and F. E. Paige, Phys. Rev. D 60, 095002 (1999) [arXiv:hep-ph/9812233]. 

[5] J. R. Ellis, A. R. Raklev and O. K. Oye, JHEP 0610, 061 (2006) [arXiv:hep-ph/0607261]. 

[6] A. Nisati, S. Petrarca and G. Salvini, Mod. Phys. Lett. A 12, 2213 (1997) [arXiv:hep- 
ph/9707376]. 

[7] G. Polesello and A. Rimoldi, ATLAS Internal Note ATL-MUON-99-006. 

[8] S. Ambrosanio, B. Mele, S. Petrarca, G. Polesello and A. Rimoldi, JHEP 0101, 014 (2001) 
[arXiv:hep-ph/0010081]. 

[9] J. Ellis, A. R. Raklev and O. K. Oye, ATLAS Note ATL-PHYS-PUB-2007-016; ATL- 
COM-PHYS-2006-093. 

[10] W. Buchmuller, K. Hamaguchi, M. Ratz and T. Yanagida, Phys. Lett. B 588, 90 (2004) 
[arXiv:hep-ph/0402179]. 



33 



[11] K. Hamaguchi, Y. Kuno, T. Nakaya and M. M. Nojiri, Phys. Rev. D 70, 115007 (2004) 
[arXiv:hep-ph /0409248] . 

[12] J. L. Feng and B. T. Smith, Phys. Rev. D 71, 015004 (2005) [Erratum-ibid. D 71, 019904 
(2005)] [arXiv:hep-ph/0409278]. 

[13] A. Rajaraman and B. T. Smith, Phys. Rev. D 76, 115004 (2007) [arXiv:0708.3100 [hep- 
ph]]. 

[14] M. Drees and X. Tata, Phys. Lett. B 252, 695 (1990). 

[15] J. L. Feng and T. Moroi, Phys. Rev. D 58, 035001 (1998) [arXiv:hep-ph/9712499]. 

[16] S. P. Martin and J. D. Wells, Phys. Rev. D 59, 035008 (1999) [arXiv:hep-ph/9805289]. 

[17] S. Dimopoulos, S. D. Thomas and J. D. Wells, Nucl. Phys. B 488, 39 (1997) [arXiv:hep- 
ph/9609434]. 

[18] S. K. Gupta, B. Mukhopadhyaya and S. K. Rai, Phys. Rev. D 75, 075007 (2007) 
[arXiv:hep-ph/0701063]. 

[19] V. D. Barger, R. W. Robinett, W. Y. Keung and R. J. N. Phillips, Phys. Lett. B 131, 
372 (1983). 

[20] H. E. Haber, arXiv:hep-ph/9405376. 

[21] B. K. Bullock, K. Hagiwara and A. D. Martin, Nucl. Phys. B 395, 499 (1993). 

[22] M. R. Buckley, H. Murayama, W. Klemm and V. Rentala, arXiv:0711.0364 [hep-ph]. 

[23] G. Corcella et al, arXiv:hep-ph/0210213. 

[24] P. Richardson, JHEP 0111, 029 (2001) [arXiv:hep-ph/0110108]. 

[25] S. Moretti, K. Odagiri, P. Richardson, M. H. Seymour and B. R. Webber, JHEP 0204, 
028 (2002) [arXiv:hep-ph/0204123]. 

[26] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12, 375 (2000) [arXiv:hep- 
ph/9903282]. 

[27] S. Jadach, Z. Was, R. Decker and J. H. Kuhn, Comput. Phys. Commun. 76, 361 (1993). 
[28] E. Richter-Was, arXiv:hep-ph/0207355. 



34 



[29] C. G. Lester and D. J. Summers, Phys. Lett. B 463, 99 (1999) [arXiv:hep-ph/9906349]. 

[30] K. Kawagoe, M. M. Nojiri and G. Polesello, Phys. Rev. D 71, 035008 (2005) [arXiv:hep- 
ph/0410160]. 

[31] H. C. Cheng, J. F. Gunion, Z. Han, G. Marandella and B. McElrath, JHEP 0712, 076 
(2007) [arXiv:0707.0030 [hep-ph]]. 

[32] M. Davids et al., CMS Note 2006/077. 

[33] K. J. F. Gaemers and G. J. Gounaris, Z. Phys. C 1, 259 (1979); K. Hagiwara, 
R. D. Peccei, D. Zeppenfeld and K. Hikasa, Nucl. Phys. B 282, 253 (1987); P. Mery, 
M. Perrottet and F. M. Renard, Z. Phys. C 36, 249 (1987); M. S. Bilenky, J. L. Kneur, 
F. M. Renard and D. Schildknecht, Nucl. Phys. B 409, 22 (1993). 

[34] P. Abreu et al. [DELPHI Collaboration], Phys. Lett. B 502 (2001) 9 [arXiv:hep- 
ex/0102041]; P. Achard et al. [L3 Collaboration], Phys. Lett. B 586, 151 (2004) 
[arXiv:hep-ex/0402036]; S. Schael et al. [ALEPH Collaboration], Phys. Lett. B 614 
(2005) 7; J. Abdallah et al. [DELPHI Collaboration], Eur. Phys. J. C 54, 345 (2008) 
[arXiv:0801.1235 [hep-ex]]. 



35 



